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Abstract 

Optical frequency conversion by four-wave mixing (Bragg scattering) in a fiber is 
considered. The evolution of this process can be modeled using the signal and 
idler amplitudes, which are complex, or Stokes-like parameters, which are real. The 
Stokes-space formalism allows one to visualize power and phase information simul- 
taneously, and produces a simple evolution equation for the Stokes parameters. 
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1. Introduction 

Parametric devices based on four-wave mixing (FWM) in fibers can amplify, frequency 
convert, pliase conjugate, regenerate and sample optical signals in communication systems 
[HIS]. The subject of this report is the nondegenerate FWM process called Bragg scattering 
(BS), in which a sideband (signal) photon and a pump photon are destroyed, and different 



sideband (idler) and pump photons are created (vr^ + vr^ 



TTp, where vr,- represents a 



photon with frequency Uj). The frequencies of the interacting waves are illustrated in Fig. 1. 
There is considerable interest in BS [3-14], because of its ability to generate an idler whose 
frequency is tunable [6], and which is not polluted by excess noise 



CO 



CO 



Figure 1: Frequency diagrams for Bragg scattering in a fiber. The long arrows denote strong pumps 
{p and q), whereas the short arrows denote weak sidebands (s and i). 

BS is described by the sideband amplitudes As and Ai. These complex amplitudes are the 
components of a Jones vector and their spatial evolution is a linear transformation in Jones 
space. Two complex variables are equivalent to four real variables. However, an inspection 
of the BS equations shows that one variable is ignorable, so BS can be described by only 
three real variables, which are the components of a Stokes-like vector. The Stokes-space 
description of BS is useful because it allows one to visualize power and phase information 
simultaneously, and because evolution in Stokes space is a simple rotation of the Stokes 
vector. 

This report is organized as follows: In Sec. 2 the coupled equations for the sideband 
(mode) amplitudes are stated and solved. In Sec. 3 these scalar equations are rewritten as 
a matrix equation for the Jones vector and their solutions are used to specify the associated 
transfer matrix. In Sec. 4 the Stokes parameters associated with the sideband amplitudes 
are defined and the (unknown) elements of a Stokes-space rotation matrix are related to the 
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(known) elements of a Jones-space transfer matrix. An evolution equation for the Stokes 
vector is also derived and solved directly. In Sec. 5, the Poisson-bracket formalism for BS is 
developed and shown to be consistent with the Stokes-vector equation. In Sec. 6, the Pauli 
spin-matrix formalism is used to relate the (unknown) elements of a Jones-space transfer 
matrix to the (known) elements of a Stokes-space rotation matrix. The mathematics of the 
bracket and spin-matrix formalisms are similar to the mathematics of angular-momentum 
operators. Hence, familiarity with the former facilitates the transition from the classical 
model of BS to the quantal model (which will be described elsewhere). Finally, in Sec. 7 the 
main results of this report are summarized. 

The evolution of a monochromatic wave with two polarization components was described, 
in Jones space and Stokes space, in [15] . The notation and results of [15] will be used without 
further comment. 

2. Coupled modes 

The initial evolution of BS (during which the pumps are not depleted) is governed by the 
Hamiltonian 

H = 5(|A|2 - + jAlA, + YAAl (1) 

together with the Hamilton equations 

d,Aj = idH/dA*. (2) 

For the case in which the pump and sideband polarizations are parallel, the wavenumber 
mismatch 6 = (/3s + — /3j — /3p)/2 + 7x(|v4pp — |y4gp)/2, where Pj = (3{ujj) are wavenumbers, 
•"yx is the Kerr coefficient and Ap and Aq are pump amplitudes, and the coupling coefficient 
7 = iSj. Other polarization configurations are discussed in [HI [12]. By combining 

Eqs. ([T]) and ([2]), one obtains the (linear) coupled-mode equations 

d^As = i6As + i-fAi, (3) 
d,Ai = i^As-iSAi. (4) 

Similar equations govern power transfer in a directional coupler [TB] and sum-frequency 
generation in a crystal [17J. Equivalent equations for the sideband powers Pj = \Aj\'^ and 
phases (pj = a.Tg{Aj) are derived in the Appendix. 
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The solutions of Eqs. ([3]) and (jlj) can be written in the input-output form 

A,{z) = /i(2;)^(0) + z/(z)A,(0), 
A,{z) = -z/*(z)^(0)+/x*(^)A,(0), 



where the transfer functions 



/i(z) = cosikz) + i5 sm{kz) /k, 
v{z) = i'ysm{kz)/k 



(5) 
(6) 

(7) 



and the BS wavenumber k = {5'^ + \'y\'^Y^'^. Notice that the transfer functions satisfy the 
auxihary equation + |z/p = 1, which is a manifestation of power conservation. The 
signal-to-idler conversion efficiency |y4j(z)p/|74s(0)p attains its maximum of |7p/A;^ when 
kz = 7t/2. 



3. Jones vector 

The scalar equations and (jlj) can be rewritten as the matrix equation 

d^A = iKA, 

where the amplitude (Jones) vector A = [As,AiY and the coefficient matrix 

5 7 



(9) 



K 



7* -5 



(10) 



Equation ([9]) has the formal solution 



Aiz) = Uiz)A{0), 



(11) 



where the transfer matrix U = e'^^. Since K is hermitian, U is unitary. Hence, it can be 
written in the Caley-Klein form 



U 



(12) 



which is consistent with Eqs. (E]) and This result is a manifestation of power conserva- 
tion. Notice that Eq. (Q is equivalent to the Hamiltonian 

H = A^KA, (13) 



together with the vector Hamilton equation 

d,A = idH/dAl (14) 



4. Stokes vector 

1/2 ' m 

Define Aj = Pj e*'^-', where Pj is a sideband power and (pj is a sideband phase. Then an 
inspection of Eq. ([T]) shows that the Hamihonian depends on the phase difference (ps — 4>i, 
but not on the total phase (j)s + 4'i- Hence, BS can be described in terms of three real variables 
(not four). For such an interaction, it is natural to introduce the Stokes-like parameters 



So 
Si 

S2 

S-s 



I A. 



I A, 



As A* + AlAi, 
ii^AgA^ AgAi 



(15) 
(16) 
(17) 
(18) 



in which the sideband amplitudes Ag and Ai play the roles of the polarization components 
X and Y [18]. 5*0 is the total sideband power, 5*1 is the power difference, and 5*2 and 5*3 



contain information about the phase difference. Notice that Sq = Sf + S2 + S^. 



By combining Eqs. ([5]) and with Eq. f[T5|) . one finds that 



^o(^) = 5o(0). 



(19) 



The total power is conserved because the Hamiltonian does not depend on the total phase, 
as discussed in the Appendix. This result implies that the Stokes vector S = {Si, S2, S3) has 
constant magnitude, so evolution in Stokes space is rotation on a sphere of radius Sq (the 
Stokes sphere). By combining Eqs. ([5]) and with Eqs. f|T6|) - f[T8|) . one finds that 



S{z) = Riz)SiO), 
where the column vector S = [Si, S2, S^Y rotation matrix 



(20) 



(21) 
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The identity + |z/p = 1 ensures that R is orthogonal. By combining Eqs. ([7]), ([8]) and 



(pTj) . one obtains the exphcit formulas 

k''Rii{d) = 6^ + {j^ + -f^)cos9, (22) 

PR^2{0) = 6-fr{l - cos 9) --fik sin 9, (23) 

k'^R^^{9) = -57,(1 - cos 0) -7^A; sin ^, (24) 

k'^R2i{9) = 5jril-cos9) +-fik sin 9, (25) 

k^R22{9) = 7' + (7- +5')cose, (26) 

k'^R2s{9) = --frliil -cos 9) + 5k sin 9, (27) 

k'^R3i{9) = -5j,{l-cos9) +-frksin9, (28) 

k'^Rs2{9) = -jrliil- cos 9) - 5k sin 9, (29) 

A;^i?33(^) = 7- +('5' + 7') cos ^, (30) 



where the distance parameter 9 = 2kz. The evolution of the Stokes vector is illustrated in 
Fig. 2. 




Figure 2: Trajectory of the Stokes vector for 6/'yr = 0.3 and 74/7^ = 0.1. The thick blue and 
red lines denote the rotation axis r and input Stokes vector 5(0), respectively. The thin red curve 
denotes the trajectory of the tip of the Stokes vector. 

In the preceding analysis, known results for the mode amplitudes were used to deduce 
results for the Stokes parameters. However, one can also determine the evolution of the 
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Stokes parameters directly. By combining Eqs. ([3]) and (jl]), one obtains the evolution 
equations 

4^0 = 0, (31) 

4^1 = -27,^2-27,53, (32) 

4^2 = 2553 + 27,^1, (33) 

4^3 = 2-frSi-25S2. (34) 

Equation fl^ is consistent with Eq. f[T^ . and Eqs. fl52]) - flM|) can be rewritten as the vector 
equation 

dj=nx S, (35) 

where the rotation vector Q = 2{—6, -7^,71). [If one were to replace 7 by 7*, the rotation 
vector would be —2(5, 7^,71) and some formulas would be simpler. However, the appearance 
of 'jAi in the equation for Ag is standard.] Define Q = 2kf, where the wavenumber k was 
defined after Eq. ([H]) and r = {—6, — 7,, 7j)/A; is a unit vector parallel to Q. Then S\\ = r{f-S) 
is parallel to the rotation vector, S± = S — f{f- S) is perpendicular to the rotation vector and 
lies in the plane defined by r and 5, and the auxiliary vector S'x = r x 5 is perpendicular to 
the rotation vector and the rs-plane. If the Stokes vector rotates about r through the angle 
6 = 2kz, 5|| remains constant and S± becomes cos6'Sj_ + sin 6' S'x. By combining these facts, 
one obtains the rotation formula 

S{9) = [cos + (1- cos ^)ff-+ sin ^fx]5(0). (36) 

Equation (136|) is consistent with Eqs. (I22|) - (l30i) . The quantity in square brackets is the 
rotation operator, written in dyadic (rather than matrix) form. 

5. Poisson-bracket formalism 

Another inspection of Eq. shows that the Hamiltonian 

H = 6S, + -frS2--f^Ss, (37) 

where S1S3 were defined in Eqs. f|T6|) - f|T8|) . This equation can be rewritten in the compact 
form H = — i1 ■ S/2, where was defined after Eq. fl35|) . Since the Hamiltonian depends 
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only on the Stokes parameters, it is natural to formulate the interaction in Stokes space. For 
any Stokes parameter Sj, the rate of change 

( dS, dA, dS, dAl\ 



By combining Eqs. ([2]) and fl38|) . one obtains the Hamilton equation 

dzS,=t{S,,H}, (39) 

where the Poisson bracket 

dP dQ dP dQ 



(40) 



dAkdAl dAldAkJ 

Notice that {Q,P} = —{P,Q}. Since H depends on Sk, to determine the consequences of 
Eq. fl39l) one must first calculate the Poisson brackets {Sj, Sk}. The results are 

{S,,Sk} = ±2tSi, (41) 

where the + (— ) sign applies if j, k and / are in positive (negative) cyclic order. By combining 
Eqs. (l39l) and pil) . one obtains the evolution equations 



dzSl = ^2S3 — ^^3'S'2, 


(42) 


dzS2 = — ^^l5'3, 


(43) 


dzS^ = ^2152 — ^2Si, 


(44) 



which are equivalent to Eq. (1351) . 5*0 is constant because {Sq, Sj} = 0. The Poisson-bracket 
formulation of classical mechanics is described in 

6. Spin-matrix formalism 

In the preceding sections the Stokes-space evolution was described in terms of Jones- 
space transfer functions [Eqs. ([7]) and ([8])], and directly in Stokes space [Eq. (l35l) ]. In this 
section the Jones-space evolution is described in terms of Stokes-space quantities. 

First, suppose that a Stokes vector is written in the polar form 

S = So (cos 6', sin 6' cos (j),sm9 sin (p) . (45) 



Then the associated Jones vector 

A = Sl'''[cos{e/2)e-'^'\ sm{e 1 2)e'^ (46) 

One can verify Eq. (H6l) by combining it with Eqs. (fT6!) - (fT8l) and comparing the results to 
Eq. (H5i) . Polar coordinates are illustrated in Fig. 3. 




Figure 3: Polar coordinates in Stokes space. The red line denotes the Stokes vector S. 

Second, suppose that a Stokes-space rotation is specified by the rotation axis r and angle 
6 (which should not be confused with the polar angle of the preceding paragraph). Some 
analysis is required to relate these quantities to their counterparts in Jones space. In Sec. 3 
it was shown that the transfer matrix U = e^^^, where K is the coefficient matrix ffTOj) . One 
can facilitate the exponentiation of iKz by defining the identity matrix ctq and the Pauli 
spin matrices 





1 




1 




-i 


CTi = 


-1 




1 




i 



These hermitian matrices have the properties = ctq and ajak = ^icri, where the + (— ) 
sign applies if j, k and I are in positive (negative) cyclic order. The latter result implies that 
(Jjak + (ykC^j = 0. Furthermore, any complex matrix M can be written in the form 

M = moCTo + m ■ a, (48) 

where the scalar and vector coefficients mo = Tr(croM)/2 and rh = Tr(aM)/2, respectively. 
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and the spin vector o = (cri, ct2, da). The coefficient matrix has the decomposition 

K = 5ai + 7^0-2 - 7icr3, (49) 

which can be rewritten in the compact form K = —kf- a, where k and r were defined after 
Eq. (!35l) . By using the aforementioned properties of the spin matrices, one finds that 

U{z) = aocos{kz) — ir ■ a sm{kz). (50) 

By combining Eqs. f|T2l) and fHSl) . one also finds that 

U = fIr<Jo + iifJ'iCri + + ^rCTs)- (51) 

It follows from Eqs. and ([501) that 

fir{z) = cos{kz), [fii{z),h'i{z),h'r{z)] = — sm{kz)f, (52) 

where [fii, z/j, z/^] is a row vector. Equations (!52l) relate the Jones-space transfer functions 
to the Stokes-space rotation axis r and the distance parameter kz. The discussion between 
Eqs. ( l35l) and ( l36l) shows that fc^; = 6'/2, where 6 is the Stokes-space rotation angle. 

The connections between the formulas that describe evolution in Jones and Stokes space 
are not accidental. It follows from Eqs. (|T6l) - f|T8|) and (H7j) that 

Sj = A^ajA. (53) 

By combining Eqs. ([9]) and ( l53l) . one finds that 

d,Sj =iA^[aj,K]A, (54) 

where the commutator 

[P, Q] = PQ- QP. (55) 
The spin matrices satisfy the commutation relations 

[aj,ak] = ±2iai, (56) 

where the + (— ) sign applies if j, k and I are in positive (negative) cyclic order. Equations 
( l56l) are the analogs of Eqs. (HTi) and the identity [a,f- a] = 2ifx a follows from them. By 
combining this identity with Eq. flMj) . one obtains the evolution equation 

dj = 2kf X S, (57) 
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which is equivalent to Eq. (!35|) . 5*0 is constant because [o"o,o"j] = 0. Other connections 
between Jones space and Stokes space are described in [15] . 

7. Summary 

Optical frequency conversion by Bragg scattering (BS) in a fiber was considered. The 
evolution of BS can be modeled using the signal and idler (sideband) amplitudes [Eqs. ([3]) 
and which are complex, or the sideband powers and phases [Eqs. (1601) - (!63l) ]. which are 
real. The amplitudes are the components of a Jones vector and their spatial evolution is a 
linear transformation in Jones space [Eq. f|TT]) ]. However, the Hamiltonian for BS [Eq. ([1]) or 
fl58|) ] does not depend on the total sideband phase. This fact has two important consequences. 
First, BS can be described by only three real variables, which are the components of a Stokes- 
like vector. The Stokes-space description of BS is useful because it allows one to visualize 
power and phase information simultaneously. Second, the total sideband power is conserved. 
In Jones space the norm of the Jones vector is constant, so the transformation is unitary [Eq. 
(|T2l) ]. whereas in Stokes space the length of the Stokes vector is constant, so its evolution is 
rotation [Eq. fl36|) ]. 

The analysis of BS proceeded in three phases. In the first phase, the coupled equations 
for the sideband (mode) amplitudes were stated [Eqs. ([3]) and (jlj)] and solved [Eqs. ([7]) and 
([H])]. These scalar equations were rewritten as a matrix equation for the Jones vector [Eq. 
Q] and their solutions were used to specify the associated transfer matrix [Eq. (fT2|) ]. The 
Stokes parameters associated with the mode amplitudes were defined and the (unknown) 
elements of a Stokes-space rotation matrix were related to the (known) elements of a Jones- 
space transfer matrix [Eq. fl2Tl) ]. In the second phase, an evolution equation for the Stokes 
vector was derived [Eq. fl35l) ]. which confirmed that evolution in Stokes space is rotation. 
This equation has a concise solution [Eq. fl36|) ]. which relates the rotation axis and angle 
to the mismatch and coupling coefficients, and the distance. The Poisson-bracket formalism 
for BS was developed [Eq. (139!) ] and shown to be consistent with the Stokes-vector equation 
[Eqs. (H2l)- (H^ ]. In the third phase, the Pauli spin-matrix formalism was used to solve 
the Jones- vector equation [Eq. (I3UI) ]. This solution is consistent with the solutions of the 
coupled-mode equations and relates the (unknown) elements of a Jones- space transfer matrix 
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to the (known) elements of a Stokes-space rotation vector [Eq. (!52l) ]. Thus, exphcit formulas 
were derived, which relate the Jones and Stokes pictures of BS. The power-phase formulation 
of the BS equations is described in the Appendix. 

The mathematics of the bracket and spin-matrix formalisms are similar to the mathe- 
matics of angular-momentum operators. Hence, familiarity with the former facilitates the 
transition from the classical model of BS to the quantal model (which will be described in a 
future report). 
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Appendix: Phase plane 

Section 2 was based on the complex formulation of the Hamilton function ([1]) and equa- 
tions ([2]). This appendix is based on the associated real formulation. Define Aj = P^^'^e^'^^ 
and 7 = |7|e*'^^. Then the initial evolution of BS is governed by the Hamiltonian 

H = 6{Ps - P,) + 2|7|(P,P,)i/2 cos(0, - 0, - 0^), (58) 

together with the Hamilton equations 

d,Pj = -dH/d(f)j, d,(f)j = dH/dPj. (59) 

The assumption that 7 is real is equivalent to the assumptions that 0^ and 0j are measured 
relative to 0^/2 and —0^/2, respectively. By combining Eqs. (|58l) and (|59l) . one obtains the 
power equations 

4Ps = 27(P,Pi)i/2sin(0,-0i), (60) 
4P. = -2^iPsP^)'/^sm{(|),-(|),) (61) 

40s = 5 + 7(^^/i's)'/'cos(0,-0,), (62) 
40, = -S + j(PjP,f/\osi(Ps-(pi). (63) 

It follows from Eqs. fIBU]) and flMl) that the total power P^ = P, + Pj is conserved (because 
H depends on 0^ — 0j). It also follows from Eqs. (!60l) - (!63l) that H is constant (because it 
does not depend explicitly on z). 

Define the power difference p = {Ps—Pi)/Pt, phase difference = (ps—'Pi and Hamiltonian 
h = H/ i^Pt)- In the notation of Sec. 4, p = Si/ Sq and is the angle between the projection 
of S on the 23-plane and the 2-axis (Fig. 3). Furthermore, let 5/7 ^ 5 and 272; z. Then 
BS is governed by the normalized Hamiltonian 

/i = 5p+ (1 -p2)i/2cos0, (64) 

together with the normalized Hamilton equations 

dzP = —dh/d(f), dz(j) = dh/dp. (65) 
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and the phase equations 



By combining Eqs. (IMI) and (165!) . one obtains the normalized power and phase equations 

d,p = (1 -p^)^/2sin<^, (66) 
40 = 5-pcos0/(l-/)^/^ (67) 

respectively. Equations ( l66l) and ( 1671) are consistent with Eqs. (I60l)-( l63!) . 

Phase diagrams associated with the Hamiltonian ( l64l) are shown in Fig. 4. The phase 
point {p, 0) moves in such a way that an observer moving with it, and looking forward, 
keeps higher energies on his left. For 6 = 0.0 the trajectories are librations. The trajectory 
that starts at the point {—tt/2, 1) is a straight line to {—tt/2, —1), followed by a phase jump 
to (7r/2, —1), followed by a straight line to (vr/2, 1), followed by a jump back to (— 7r/2, 1). 
For 6 = 0.3 some trajectories are librations, whereas others are rotations. The trajectory 
that starts at (— 7r/2, 1) is a curve to the point {tt/2, 1), followed by a phase jump back to 
(— 7r/2, 1), and the trajectory that starts at (— 7r/2, — 1) is a curve to the point (7r/2, — 1), 
followed by a jump back to (— 7r/2, —1). 




-3-2-10123 -3-2-10123 
Phase Phase 



Figure 4: Phase diagram for (a) S/j = 0.0 and (b) S/j = 0.3. Lighter regions correspond to higher 
energies, whereas darker regions correspond to lower energies. The energy contours (solid curves) 
are the trajectories of the phase point (p, 0). 

These phase diagrams describe the evolution of BS qualitatively. To describe the evo- 
lution quantitatively, one must solve Eqs. (jHSD and (1^ . First, suppose that 6 = 0, which 
corresponds to (left-handed) rotation about the 2-axis in Stokes space. Then 

(4p)' = (l-p')sin2 0, 

= l-h^-p"^. (68) 
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It follows from Eq. (!68l) that 

p{z) = {l-h^y/^cos{z + zo), (69) 

where h is constant and the distance parameter zq = cos^^[p(0)/ (1 — h'^Y^'^]. For the case in 
which 5 = and ji = 0, Eqs. (1201) . ( l22l) and can be rewritten in the normalized form 

= Si(0) cosz — 53(0) sin^;. (70) 

By using the identity sl + sf = 1 — s^, one finds that Si(0) = (1 — s'^Y^'^ cos and 53(0) = 
(1 — s^)^''^ sin Zq, where S2 = h and zq = tan~"'^[s3(0)/si(0)] is the angle between the projection 
of 5'(0) on the 31-plane and the 1-axis. Hence, 

s^{z) = {l-slY/'cos{z + zo). (71) 

Equation (!7T]) is equivalent to Eq. (169!) . 

Second, suppose that 6^0, which corresponds to (left-handed) rotation about the axis 
(5, 1, 0)/k, where k = {1 + 6^/^ Then 

(4p)' = l-p'-{h-6pY, 

= k'[pl-ip-pan (72) 

where the power average Pa = Sh/k"^ and power difference pd = {k"^ — h?)^/"^ /k"^. In the 
notation of Sec. 4, pa = ri(r ■ S) and pd = |Sj_(0)|/A;. It follows from Eq. (1721) that 

p{z) = pa+PdCos[k{z + Zo)], (73) 

where kzo = cos'^{[p(0) —Pa]/Pd}- For the case in which 7^ = 0, Eqs. (l20l) and (I22l) - (l2^ 
can be rewritten in the normalized form 

k'^siiz) = si{0)[5^ + cos{kz)] + S2{0)5[l-cos{kz)]-S3{0)ksm{kz), 

= 5[SsiiO) + S2iO)] + [siiO)-Ss2iO)]cos{kz)-ks3iO)sm{kz), (74) 

where 5si(0) + ^2(0) = h. By using the identity + + S3 = 1, one finds that [si(0) — 
5s2i0)]^ + [kss{0)]^ = k^-h"^. Hence, Eq. ([71]) can be rewritten in the form of Eq. ([73]), with 
kzo = tan~^{A;s3(0)/[si(0) — 5s2(0)]}. The nonlinear evolution of BS is described in [Slin]- 
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